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l.  Introduction 


Fractional  factorial  designs  are  among  the  most  utilized  statistical  designs  by  prac- 
tioners.  In  a  series  of  three  papers  Rao  (1946,  1947,  1949)  identified  special  types  of 
fractions  with  a  great  deal  of  symmetry  and  many  desirable  statistical  properties.  These 
special  fractions  are  now  known  as  orthogonal  arrays.  Numerous  beautiful  results  have 
been  obtained  by  both  statisticians  and  mathematicians  on  this  topic.  These  results  are 
published  in  a  wide  variety  of  journals  and  presented  in  many  different  styles.  In  their 
forthcoming  research  monograph  Hedayat  and  Stufken  (1986)  have  presented  the  various 
results  on  this  fascinating  subject  in  a  unified  and  comprehensive  way. 

r 

This  paper  is  divided  into  four  sections.  In  section  2.we  will  review^basic  definitions 
and  terminology  of  the  subject.  /In  section  3-wewill  present  an  efficient  way  of  preparing 
design  and  information  matrices  associated  with  orthogonal  arrays  under  orthogonal  poly¬ 
nomial  models.  Section  4  contains  a  new  result  and  a  new  concept.  The  new  result  states 
that  with  orthogonal  arrays  of  strength  t  we  can  orthogonally  estimate  other  parametric 
vectors  besides  the  one  which  is  advocated  in  the  literature.  We  have  also  identified  some 
types  of  orthogonal  arrays  which  are  somewhere  between  orthogonal  arrays  of  strength  t 

”  ‘  t  '  ?  a  *  c 

and  t  +  1.  We  caff  them  flexible  orthogonal  arrays  of  strength  t.  Practical  applications  of 
such  arrays  are  also  pointed  out.  These  results  are  very  useful  in  preparing  line  graphs  of 
Taguchi  for  orthogonal  arrays. 


2.  Preliminaries 


We  begin  by  giving  some  of  the  basic  definitions  and  terminology.  Let  5  be  a  set  of 
s  symbols,  coded  by  0, 1, . . .  ,s  —  1.  A.  k  x  N  array  A  with  entries  from  S  is  called  an 
orthogonal  array,  denoted  by  OA(N,k,s,t),  if  each  t  x  N  subarray  of  A  has  the  property 
that  every  possible  (  x  1  vector  with  entries  from  S  appears  equally  often  (say  A  times)  in 
the  columns  of  the  chosen  subarray.  The  integer  A  is  called  the  index  of  the  array,  while 
N,k,s  and  t  are  said  to  be  the  parameters  of  the  array.  If  A  =  1  the  array  is  said  to  be  of 
index  unity.  The  relation  N  =  As*  is  an  immediate  consequence  of  the  definition. 

The  reader  may  verify  that  the  following  array  is  an  OA (8,4 ,2,3),  an  orthogonal  array 
of  index  unity. 
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01111000 
10  110  10  0 
110  10  0  10 
1  1  1  0  0  0  0  1 

The  terminology  that  is  used  in  the  literature  for  the  parameters  of  an  orthogonal  array 
is  as  follows:  N  is  called  the  size  of  the  array  or  the  number  of  runs,  assemmblies  or 
treatment  combinations;  k  is  called  the  number  of  constraints  or  factors  in  the  arrays;  s  is 
called  the  number  of  levels  or  symbols;  t  is  called  the  strength  of  the  array. 

3.  ORTHOGONAL  ARRAYS  UNDER  ORTHOGONAL  POLYNOMIAL  MODELS 

As  a  fractional  factorial  designs  orthogonal  arrays  are  highly  efficient  if  orthogonal 
polynomial  models  are  postulated  for  the  responses  under  study.  It  was  indeed  in  this 
context  that  Rao  introduced  the  concept  of  orthogonal  arrays  in  the  literature.  The  purpose 
of  this  section  is  two-fold.  To  present  the  orthogonal  polynomial  models  associated  with 
orthogonal  arrays  so  that  our  result  in  section  4  will  be  easily  understood  by  the  readers  of 
this  paper.  We  also  present  an  efficient  way,  recently  discovered  by  Hedayat  and  Stufken 
(1986),  for  generating  the  design  and  information  matrices  associated  with  orthogonal 
arrays.  This  new  method  is  very  useful  in  analyzing  data  obtained  under  orthogonal 
arrays  and  in  studying  the  optimality  of  such  arrays. 

In  the  remainder  of  this  section  we  will  have  to  assume  that  the  reader  has  some 
familiarity  with  the  standard  concepts  in  factorial  experiments  and  the  theory  of  ANOVA. 
For  a  quick  review  of  these  concepts  a  useful  reference  is  Raktoe,  Hedayat  and  Federer 
(1981).  Assume  that  our  interest  is  in  an  experiment  based  on  k  controllable  factors, 
each  of  which  can  be  set  at  s  different  levels.  The  levels  are  coded  by  0,1,..., s  -  1. 
A  treatment  combination  is  an  experimental  condition  to  which  each  factor  contributes 
a  level.  For  example  (t'i, . . . ,  ij, . . .  ,**)  is  a  treatment  combination  with  the  jth  factor 
at  level  ij.  Thus  there  are  sk  distinct  treatment  combinations.  These  can  be  exhibited 
column  wise  in  a  Ik  x  s*  array,  which  we  will  denote  by  p  and  refer  to  as  a  minimal  complete 
factorial  based  on  k  factors  at  s  levels  each. 

Let  g  be  a  treatment  combination  and  yg  the  observed  response  under  the  experimental 
condition  g.  We  will  assume  that  yg  is  a  random  variable, 
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y*  =  f'W  +  £g, 

where  f'(g )  is  a  row  vector  of  real  known  functions  of  g  and  9  is  a  column  vector  of 
unknown  parameters.  The  random  error  component  eg  will  be  assumed  to  have  mean  zero 
and  unknown  variance  a2. 

The  structure  of  f'(g)  is  an  important  consideration  in  practice.  If  all  factors  are 
quantitative  factors  a  popular  way  of  modeling  yg  is  to  structure  f'(g)  via  the  orthogonal 
polynomial  model  (see  for  example  Chapter  4  of  Raktoe,  Hedayat  and  Federer  (1981)). 
We  will  use  this  model  from  here  on.  If  we  perform  the  experiment  using  all  sk  treatment 
combinations,  the  resulting  vector  of  observations,  Yp,  can  be  expressed  as 


(3.1)  Yp  —  Xp0p  +  €p, 

where  Xf  is  the  design  matrix  and  /3P  is  the  vector  of  general  mean,  main  effects  and 
interactions.  The  entries  in  0P,  following  the  common  notation,  are 

4% 4%  •••  4%  »  *;  €  {0,1 . a  -  1}. 

We  will  now  present  an  easy  way  to  find  the  entries  of  Xp  under  the  orthogonal  polynomial 
model.  This  model  is  especially  useful  if  we  axe  only  interested  in  a  fraction  of  the  minimal 
complete  factorial  design.  Let  Xo,Xi,...,xa_i  be  the  actual  levels  of  the  first  factor.  Form 
the  following  matrix 


-1 

xo 

x0 

T*-1 

—  xo 

1 

X\ 

X? 

X*-1 

Xj 

Xi  = 

• 

; 

.1 

X,-l 

xs-l 

H 

tm  fe 

1  1  * 

M  M 

Orthogonalize  the  columns  of  X\  from  left  to  right  and  call  the  resulting  matrix  M\. 
Similarly  obtain  M2, . . .  ,M*_  1  corresponding  to  the  other  k  —  1  factors.  Denote  the  entry 
in  the  ith  row  and  jth  column  of  M<  by  njj\  0  <  j  <  a  -  1,  1  <  l  <  k.  Then  the 
entry  in  Xp  corresponding  to  the  treatment  combination  (x,, ,  xtJ , . . . ,  x,k )  and  parameter 
tfx  ^  4>3k  is  8iven  by 


(3.2) 


n(*>  -n(2).  n(k) 

nUJl  nt2}2  ‘  '  ‘  nik]k  ' 


If  the  5  levels  of  the  factors  are  equally  spaced,  the  orthogonalized  matrices  M  can  be 
obtained  from  the  available  literature.  Table  4.1  in  Raktoe,  Hedayat  and  Federer  (1981) 
lists  such  matrices  for  2  <  s  <  7,  which  is  sufficient  for  most  practical  purposes.  The  same 
table  can  be  used  for  s  =  8  or  9  if  a  polynomial  model  of  degree  less  than  or  equal  to  5  is 
fitted  for  each  factor.  We  demonstrate  some  of  the  above  ideas  in  the  following  example. 

Example  3.1.  Let  s  =  Z,k  =  2.  Under  equally  spaced  levels,  coded  by  0, 1,2,  we  have 


P  = 


000111222 
0  1  2  0  1  2  0  1  2 


Mi  =  Mi  = 


0  12 

/i  -i  i' 
1  0  -2 
\1  1  1 J 


Pp  —  ^1$2»  ^1^2)  ^1^2>  4>\4>2i  4>\4,2i  ^l^)' 


Thus  for  example,  the  entry  in  Xp  corresponding  to  treatement  combination  (1,2)  and 
parameter  is,  from  (3.2)  equal  to 


»12»20  =  (~2)  •  l1)  =  “2- 

Completing  the  example  we  obtain 
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It  is  often  impractical  or  impossible  due  to  other  restrictions  (money,  time,  etc.)  to 
use  all  possible  treatment  combinations  in  an  experiment.  Clearly  we  would  like  to  select 
a  fraction  of  the  sk  treatment  combinations,  where  we  will  allow  that  some  of  the  selected 
combinations  axe  used  more  than  once  in  the  experiment.  If  A  is  an  OA(N,k,s,t)  based  on 
5  =  {0, 1, . . . ,  s  —  1},  then  the  columns  of  A  form  such  a  fraction  of  the  minimal  complete 


factorial. 


Under  the  orthogonal  polynomial  model  (3.1),  possibly  neglecting  some  of  the  param¬ 


eters  in  0P,  we  can  write 


YA  =  XA0A  +  eA. 


Here  0A  is  the  vector  of  parameters  obtained  from  0P  by  deleting  those  parameters  that 
can  be  neglected,  while  XA  is  the  design  matrix  obtained  as  in  (3.2),  but  now  only  using 
the  treatment  combinations  in  A  and  the  vector  of  parameters  0A. 


Example  4.1.  Let  A  be  the  following  OA(8,4,2,3): 


0  11110  0  0 
10  110  10  0 
110  10  0  10 
1  1  1  0  0  0  0  1 


Let  0A  consists  of  the  general  mean,  the  main  effects  and  the  two  factor  interactions 
between  factors  1  and  2,  1  and  3,  and  1  and  4.  In  our  notation  this  means 


0A  —  (^102^3^4)  ^1^2^3^4>  $1$2^3^4»  ^1  <$2<?3<^4> 


Since  M\  =  M2  =  M3  =  M4  = 


we  obtain  by  (3.2): 
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1-1  1  1  1-1  -1  -1 

11-111-111 

111-111-11 
1111-111-1 
1  1-1-1  -1  -1  -1  -1 

1-1  1  -1  -1  -1  1  1 

1-1-1  1-1  1-1  1 
.1  -1-1-1  1  1  1  -1. 


0A  + 


It  is  interesting  to  observe  that  in  the  above  example  X'AX\  =  81,  implying  that 
the  general  mean,  the  main  effects  and  the  3  two  factor  interactions  can  be  estimated 
orthogonally,  assuming  that  other  parameters  can  be  neglected.  This  property  is  not  based 
on  coincidence.  We  will  soon  see  a  result  that  explains  this  property.  In  the  meantine  the 
reader  is  invited  to  verify  that  the  same  property  holds  if  we  replace  the  interaction  between 
factors  1  and  4  in  /?x  by  the  interaction  between  factors  2  and  3. 

We  will  proceed  by  giving  some  of  the  statistical  properties  of  orthogonal  arrays 
as  fractional  factorial  designs  under  the  orthogonal  polynomial  model,  but  for  a  concise 
formulation  we  need  one  more  concept.  If  r  is  even,  then  a  fractional  factorial  design  is  said 
to  be  of  resolution  r  if  all  factorial  effects  involving  (r  —  2) /2  or  fewer  factors  are  estimable, 
assuming  that  all  factorial  effects  involving  (r  +  2)  /2  or  more  factors  can  be  neglected.  If  r 
is  odd,  the  design  is  said  to  be  of  resolution  r  if  all  effects  involving  (r  — 1)/2  or  fewer  factors 
are  estimable,  assuming  that  the  remaining  effects  can  be  neglected.  The  usual  convention 
about  the  general  mean  is  that  for  odd  r  it  is  among  the  parameters  to  be  estimated,  while 
for  even  r  it  is  not  among  those  to  be  estimated,  nor  among  those  that  can  be  neglected. 
So  for  example  a  resolution  III  design  is  one  for  which  the  general  mean  and  main  effects 
can  be  estimated,  assuming  that  all  other  effects  can  be  neglcted.  A  resolution  IV  design 
allows  us  to  estimate  all  main  effects  under  the  presence  of  the  general  mean  and  all  two 
factor  interactions,  while  all  other  effects  are  neglected. 

The  main  reason  that  orthogonal  arrays  axe  interesting  in  the  above  described  statis¬ 
tical  context  is  due  to  the  following  well  known  result. 

THEOREM  4.1.  An  orthogonal  array  of  strength  t  is  a  design  of  resolution  t- 1-1.  Moreover 
the  concerned  effects  are  orthogonally  estimable,  while  the  general  mean  can  always  be 
included  in  the  effects  to  be  estimated,  both  for  even  and  odd  t. 
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A  proof  of  this  result  can,  for  example,  be  found  in  Raktoe,  Hedayat  and  Federer 
(1981)  (Theorem  13.1).  This  is  a  very  nice  and  powerful  result.  It  tells  us  for  example 
that  in  an  OA(8,4,2,3)  the  general  mean  and  main  effects  can  be  estimated  orthogonally, 
assuming  that  three  and  four  factor  interactions  are  absent.  However  a  property  as  in 
Example  3.2  cannot  be  explained  by  this  result. 

For  this  purpose  we  formulate  here  a  more  general  result.  Let  x  be  a  fc-dimensional 
(0,  l)-vector,  i.e.,  x'  =  (xi,  22, . . . ,  2*)  with  2,  6  {0,1}.  With  factorial  effects  corre¬ 
sponding  to  x  we  will  mean  those  effects  <f>\1  ^  . . .  <f>l£  for  which  ty  =  0  if  and  only  if 
Xj  =  0, 1  <  j  <  k. 

Now  assume  that  B  and  C  are  two  disjoint  sets  of  k-dimensional  (0,  l)-vectors  (C  may 
be  the  empty  set)  such  that 

1.  If  x, y  €.  B,  then  |{*  6  {1, . . . , k}  :  x,-  =  1  or  yt-  =  l}|  <  t. 

2.  If  x  6  B,y  €  C,  then  |{*  €  {1, . . . ,  k}  :  Xi  =  1  or  y,  =  l}|  <  t. 


Then  we  have  the  following  theorem. 


THEOREM  4.2.  If  A  is  an  orthogonal  array  of  strength  t,  all  effects  corresponding  to 
vectors  in  B  can  be  estimated  orthogonally  assuming  that  effects  corresponding  to  vectors 
in  C  are  the  only  other  ones  that  are  not  neglected.  The  zero  vector,  corresponding  to  the 
general  mean,  can  always  be  chosen  in  B. 


It  is  easy  to  see  that  this  is  a  generalization  of  Theorem  4.1.  Moreover  if  we  choose  B  = 

{(0,0, 0,0)',  (1,0,0,0V,  (0, 1,0,0)',  (0,0, 1,0)',  (0,0,0, 1)',  (1, 1,0,0)',  (1,0, 1,0)',  (1,0,0, 1)'} 

and  C  =  <j>  we  obtain  an  explanation  for  Example  3.2. 


More  information  on  the  structure  of  the  orthogonal  array  A  can  even  lead  to  more 


general  conclusions.  Instead  of  a  general  formulation  we  just  illustrate  this  by  an  example. 


Example  4.3.  Let  A  be  the  following  orthogonal  array 


Factor 

1 

2 

3 

4 

5 


0  0  0  0  1  1  1  1 
0  0  1  1  0  0  1  1 
0  10  10  10  1 
0  10  110  10 
0  110  10  0  1 


r  •  *  j 

& 
L-'-  - 


This  is  an  OA(8,5,2,2).  We  claim  that  the  following  8  effects  can  be  estimated  orthogonally, 
assuming  that  no  other  effects  are  present:  the  general  mean,  the  main  effects,  and  the  two 
factor  interactions  between  the  factors  1  and  2  and  1  and  5.  We  can  not  conclude  this  from 
our  previous  result;  indeed  we  do  not  claim  that  our  statement  is  valid  for  any  OA(8,5,2,2). 
only  for  this  particular  one.  An  explanation  for  the  vailidity  of  our  statement  is  as  follows: 
select  any  two  effects  to  be  estimated  and  look  at  the  rows  in  A  corresponding  to  the  factors 
that  are  involved.  These  rows  form  each  time  one  or  more  copies  of  a  minimal  complete 
factorial.  For  example,  if  we  select  the  main  effect  of  2  and  the  interaction  between  1 
and  5,  we  see  that  the  rows  corresponding  to  the  factors  1,  2  and  5  form  one  copy  of  a 
minimal  complete  factorial.  Similar  for  other  effects.  Due  to  this  property  the  validity  of 
our  statements  follows.  We  call  orthogonal  arrays  with  such  additional  features  as  flexible 
orthogonal  arrays.  The  study  of  the  existence  and  construction  of  flexible  orthogonal  arrays 
is  currently  under  investigation. 

Closing  Remarks.  Our  Theorem  4.2  and  the  notion  of  flexible  orthogonal  arrays  are 
useful  in  preparing  and  efficiently  cataloging  line  graphs  associated  with  orthogonal  arrays. 
The  usefulness  of  such  graphs  for  practical  applications  of  orthogonal  arrays  is  nicely 
demonstrated  by  Taguchi  and  Wu  (1979). 
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